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ABSTRACT 



The effects of Gaussian-shaped IF and video filters on the 
performance of a square-law detector employing post-detection 
integration are analyzed. The number of additional noise-only 
samples that are integrated due to a finite video bandwidth is 
determined. Emerson's method is used to obtain an expression 
for the cumulants of the output probability density function. 
These cumulants are used in Edgeworth's asymptotic series 
expansion of the density functions. By integrating these 
density functions, the receiver operating characteristics are 
determined for various ratios of IF bandwidth to video 
bandwidth. A collapsing ratio which takes into account the 
Gaussian shape of the filters is formulated and compared 
against Barton's approximate formula for the collapsing ratio. 
For typical video bandwidths, Barton's approximate formula is 
found to overestimate the collapsing loss by an amount less 
than 0.5 dB. 
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I . INTRODUCTION 

This report is concerned with the effects of insufficient 
video bandwidth on the detection performance of a square law 
detector with Gaussian-shaped intermediate frequency (IF) and 
video filters in pulsed radar applications. 

The square-law detector is assumed to have an arbitrary 
wide low-pass zonal filter [Ref. 1] so that its instantaneous 
output is just the square of the envelope of the input signal. 
This assumption is in line with Marcum's definition of the 
detector as "any device whose instantaneous output is a 
function of the envelope of the input wave only" [Ref. 2]. 

In most practical radars with superheterodyne receivers, 
the second detector is usually followed by one or more stages 
of video amplification prior to post-detection integration and 
threshold detection. The bandwidth of the video amplifier 
must be sufficiently wide to avoid distortion of the input 
video signal from the detector. It has been pointed out that 
the video bandwidth B v must be greater than one-half the IF 
bandwidth £ lf to be "sufficiently wide" [Ref. 3], In 
practice, the video bandwidth may be less than half the IF 
bandwidth. 

An insufficient video bandwidth will have two effects on 
the detection statistics of the radar receiver. Firstly, it 
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will affect the threshold setting for a specified probability 
of false alarm, P fa . Secondly, it will adversely affect the 
input signal-to-noise ratio required to attain a desired 
probability of detection, P d , for a specified probability of 
false alarm. An equivalent statement to the second effect is 
that there will be a degradation in probability of detection 
vis-a-vis the probability of detection attainable with an 
ideal video bandwidth for the same signal-to-noise ratio and 
probability of false alarm. 

The focus of this report will be on the degradation in 
performance associated with the second effect. Mathematically, 
Marcum [Ref. 2] has treated this degradation in performance by 
modelling it as a consequence of integrating M additional 
noise-only samples along with the desired N signal-plus-noise 
samples 1 . The degradation in performance is quantified by 
introducing a "collapsing loss", L c , which has been defined in 
dB as the incremental increase in signal-to-noise ratio with 
respect to the signal-to-noise ratio required in the ideal 
case of integrating only signal-plus-noise samples for the 



'Marcum's paper first introduced the general concept of a 
collapsing loss resulting from the integration of a greater 
number of noise-only samples than signal-plus-noise samples. 
Insufficient video bandwidth is treated, among others, as a 
specific cause of collapsing loss. 
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same probability of detection. Mathematically, this is 

expressed as: [Ref. 4] 

L c dB = X N+M dB - X N dB ( 1 ) 



where X N+M is the signal-to-noise ratio required when N signal- 
plus-noise samples plus M extra noise samples are integrated, 
and X N is the signal-to-noise ratio required when only N 
signal-plus-noise samples are integrated. 

The signal-plus-noise samples and noise-only samples are 
assumed to be independent and respectively identically 
distributed. 

Mathematically, Marcum showed that for a non-fluctuating 
target, the effect on the detection statistics of integrating 
M additional noise samples can be taken into account by multi- 
plying the joint characteristic function of the N signal-plus- 
noise samples by the joint characteristic function of the M 
noise-only samples. Defining the collapsing ratio, 



Mj_N 

N 



( 2 ) 



Marcum showed that for the square law detector 



('N-fM ^ ^ ) 



exp [ - Np ( — ) ( ) ] 

P i+jS 

(l +j£)" p 



(3) 
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where C N+H ( £ ) is the characteristic function for the joint 
density function of N signal-plus-noise samples plus M 
additional noise samples. 

When p is equal to one (i.e., M is equal to zero), Eqn 
(3) reduces to the equation of the characteristic function 
with no collapsing loss, i.e., 



C N ( 5 ) 



exp [ - NX ( ) ] 



(1 +j5) " 



(4) 



Comparing Eqn (3) with Eqn (4), it is seen that the 
statistical results for p = 1 can be easily extended for any 
p. In particular, receiver operating characteristics (ROC) 
such as Meyer plots 2 for p = 1 are available in the 

literature and can be used directly to determine the signal- 
to-noise required by using pN instead of N as the number of 

x 

pulses integrated and with — as a parameter of the ROC 
instead of X. [Ref. 4] 



2 These are plots of P D as a function of the number of 
pulses noncoherently integrated for a given Fehlner's false 
alarm number and with the input signal-to-noise ratio for a 
single pulse as a running parameter. 
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Hence, the quantitative effects on the detection 
statistics of insufficient video bandwidth can be determined 
once the collapsing ratio, p, is known. The problem is in 
relating M (and hence, p) to the video bandwidth. What is 
the effective number of additional noise-only samples that are 
integrated along with the intended number of signal-plus-noise 
samples given a finite video bandwidth? 

A common practice is to use the following semi-empirical 
formula given by Barton [Ref. 5] to estimate the collapsing 
ratio : 



Pb 



Bit * 2B V 
2 B„ 



(5) 



where the IF noise equivalent bandwidth B if is assumed to be 
matched to the pulse width such that B if x = 1. The video 
filter is modeled as an ideal low-pass filter with a flat 
frequency response from -B v to +B V and zero elsewhere. The 
frequency response characteristics of the actual filter is not 
taken into account apart from computing its noise-equivalent 
bandwidth, B v . 

The purpose of this paper is to investigate the effects of 
Gaussian-shaped IF and video filters on the detection 
performance of the square-law detector. The primary 
consideration in the choice of Gaussian-shaped filters for 
this investigation is the mathematical convenience of the 
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Gaussian shape which yields analytical results in closed form. 

Although a Gaussian filter is physically not realizable, 
it has been found to be a good approximation in characterizing 
the frequency response of multistage RC amplifiers [Ref. 6], 
Hence, the results obtained will have direct application for 
square-law detectors with multistage pre-detection and post- 
detection RC filters, besides providing a relative measure of 
the accuracy of Barton's formula for the collapsing ratio. 

In the next section, the effects of video bandwidth on the 
detector's output probability density functions and operating 
characteristics will be analyzed. Specifically, a formula 
will be developed to determine the additional noise-only 
samples integrated due to a finite video bandwidth. This 
leads to the formulation of the collapsing ratio which will be 
discussed and compared against Barton's formula in Section 
III. Section IV summarizes the methods and findings of this 
paper . 
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II. DEGRADATION IN PERFORMANCE 



A. DESCRIPTION OF THE SYSTEM 

The system under consideration is illustrated in Figure 1. 
It consists of a Gaussian IF filter, a square-law detector, a 
video amplifier with a Gaussian frequency response, and a 
linear post-detection integrator. For convenience, it shall 
be referred to as the Gaussian square-law system. 




Figure 1. Block diagram of the Gaussian square-law system. 



The IF filter has a Gaussian passband characteristic with 
a center frequency of f Q and a root-mean-square bandwidth of 
P . The video amplifier also has a Gaussian passband but with 
a root-mean-square bandwidth of v centered at zero frequency. 
Since the gain of the video amplifier will have no impact on 
detection performance, it is assumed to have unity gain for 
mathematical convenience. Correspondingly, the impulse 
response functions of the IF filter and video amplifier are 
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respectively defined by: 




( 6 ) 



and 




(7) 



The signal to be detected is a pulse-modulated sinusoidal 
waveform corresponding to the echo from a target at a 
particular range from the radar. 

The input noise to the system is assumed to be additive 
white Gaussian noise (AWGN) with a two-sided power spectral 
density equal to x\/2 . The noise-equivalent bandwidth of the 
IF filter, £ if , is matched to the pulse width so that noise at 
the input to the detector can be assumed to be completely 
uncorrelated from pulse to pulse [Ref. 2]. 

B. ADDITIONAL NOISE-ONLY SAMPLES INTEGRATED 

For the hypothetical case of infinite video bandwidth, the 
video range bins will have the same resolution as the IF range 
bins, which is of the order of one pulse width. However, for 
most practical video amplifiers with finite bandwidths, the 
resolution of the video range bins will be poorer than those 
of the IF range bins. In such cases, each sampled output of 
the video amplifier will be shown in the subsequent paragraphs 



to be equivalent to the result of averaging N v independent 
samples of the detector output, where N v is a number greater 
than one. Since the detector follows a square law, each of 
the N v independent samples is equal to the square of the 
voltage envelope in the corresponding IF range bins. 

Suppose a target is present at a particular IF range bin. 
One of the N v samples will consist of signal-plus-noise. The 
remaining N v - 1 samples will consist of noise only. Hence, N v - 1 
additional noise-only samples will be integrated for every 
signal-plus-noise sample integrated. The added noise results 
in a degradation in performance that is quantified by a 
collapsing loss. 

The number of noise-only samples integrated along with the 
signal-plus-noise sample can be determined by recognizing that 
the video amplifier with a finite bandwidth behaves as an 
analog filter integrator and paralleling it with a digital 
integrator that averages N v independent samples of the 
detector output. This is illustrated in Figure 2. 

The digital integrator as shown in Figure 2(a) has the 
well-known effect of decreasing the variance of the input 
samples, g 3 , by a factor of N vl i.e., 



Var { y } 



Var Ig^-) 
K 



( 8 ) 



where y is the output of the digital integrator. 
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Square-Law 


0/ 


' £ 


Detector 




N v q i 

1=1 



(a) Digital Integrator: Var(y) = Var(q)/N v 




(b) Video filter: Equivalent N v = Var{q(t))/Var{y(t)) 



Figure 2 . Equivalent number of independent samples averaged 
by the video filter. 



Hence, the equivalent number of samples, N v , averaged by 
the video filter can be determined by evaluating: 



Var igjt ) ) 
Var {y ( t) ) 



(9) 



where 

q(t) = output of the square-law detector 
y(t) = output of the video filter (amplifier). 

The cumulants of the probability density function (pdf) f y 
at the output of the Gaussian video filter can be determined 
by the method of Emerson [Ref. 7]. When the input to the 
square law detector is a CW sinusoidal signal contaminated by 
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AWGN, the cumulants are given by: 



K 



(2a 2 ) 1 (r- 1) ! 



x 



i 



l 



( 10 ) 



[ ( 1 + 2y 2 ) 2 + 1] r - [ (l + 2y 2 ) 2 - 1] r 



i l 




_i 

4 



where 

K y/r = rth cumulant of the pdf , f y 

a 2 = average noise power at the input to the detector 
X ' = signal-to-noise ratio at the input to the detector 




The mathematical derivation of Eqn (10) is shown in 
Appendix A. For the purpose of this Section, it is sufficient 
to point out that the derivation takes into account the shape 
of the Gaussian filter and that the cumulants are in general 
related to the characteristic function of the associated 
probability density function by the following definition 
[Ref. 8]: 




( 11 ) 
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The cumulants are also related to the central moments of 



the density function. 3 In particular, the variance or second 
central moment of the density function is equal to the second 
cumulant of the density function. 

Hence, for the case of a CW sinusoidal signal, 



Variy(C)} = K. 



y. 2 



V 1 +2y 2 



1 + 2X 



/ 



l + 2y 2 1 



N l' + Y : 



( 12 ) 



To use Eqn (12) for evaluating N v in Eqn (9), X' must be 
the equivalent CW signal-to-noise ratio that will produce the 
same variance resulting from the averaging of one signal-plus- 
noise sample and N v -1 noise-only samples. 

Referring to Figure 2(a), let X be the input signal-to- 
noise ratio of the jth range bin containing the target. Then, 



y 



-k* 



— £ 

N, 



<2i 



i-i 



(13) 



and 



Var iy) 



1 i \ N v - 1 . s 

— Var (g,) + - -- - - Var {g^ , 

nI nI 



1 * 0 (14) 



3 A full treatment of the relationship between the 
cumulants and the moments of the probability density function 
can be found in Ref. 8. 
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where g* are independent and identically distributed noise- 
only samples for i = j. 

Now, g A is the output of the square-law detector prior to 
any video filtering. For analysis, this is theoretically 
equal to the output y when the video filter has an infinite 
bandwidth. Hence, 



Var {gj} = K s 



y ' 2 'y=0,X'~X 



= o 2 (1 + 2X) 



and 



V3r{q ‘ } - ^iy-O.X'-O 



(15) 



(16) 



- o 



2 



Substituting the above results into Eqn (14) and simplify- 
ing yields the following expression for the variance of y: 



Var (y) 



2o ; 

nI 



X + 



q 

N„ 



(17) 



Let X be the equivalent signal-noise-ratio in each of the 
N v range bins that will produce the same variance resulting 
from averaging the single signal-plus-noise sample with N v 
noise-only samples. 
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Then, using Eqn (12) with y = 0 and X' ~ X, it follows that 



Var iqj) = o 2 ( 1+2*) 

and 

r i i 

Var = V*r U- £ q, \ 

iy v i-i 

- — (1 + 2 *) 

N v 



(18) 



(19) 



By equating Eqn (17) and Eqn (19), it can be easily shown 



that 



X = 



_X_ 

N„ 



( 20 ) 



Hence, for an arbitrary y, it follows from Eqn (12) and 
Eqn (20) that 



Variy(t)) = K 



•y - 2 I x'-X 



( 21 ) 



Vl + 2y 2 



1 + 2 



X 

N v \ 



1 + 2y 2 
1 + Y 2 
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The equivalent number of independent samples W v averaged 
by the video filter can now be determined by using Eqn (9), 
i . e . , 



Var [q( t) ) 

Var {y ( t) } ( 22 ) 



with 



1 + 2y 2 
1 + Y 2 

and 

Var {g( t) ) = o 2 (1 + — ) 

N„ 



Var {y ( t) } 



»/i + 2 v 2 



_X 



(23) 



(24) 



Substituting Eqn (23) and Eqn (24) into Eqn (22) results 
in the following quadratic expression in N v : 



[ 



Ny + 2Z, 



1 + 2Y_ 
\ 1 + Y 2 



1 



- v^l + 2 Y 2 j N v - 2X\Jl + 2y 2 =0 (25) 



Hence, the equivalent number of independent samples 
averaged by the Gaussian-shaped video filter is given by: 



AT = ~ b t 

2 



(26) 
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where 



b 



and 



2X 



N 



1 + 2y z 
1 +y 2 



st 1+Y 2 



c - -2X v/l + 2y 2 



(27) 



(28) 



Solving Eqn (26) for y = 0 shows that W v is equal to one 
when the video bandwidth is infinitely large. This result is 
consistent with the earlier observation that the video range 
bins will have the same resolution as the IF range bins when 
the video bandwidth is infinite. 



C. RECEIVER OPERATING CHARACTERISTICS 

The collapsing loss caused by insufficient video bandwidth 
can be observed by examining the effect of integrating the 
additional noise-only samples on the output probability 
density function and operating characteristics of the Gaussian 
square law system. 

The number of signal-plus-noise samples N integrated by 
the post-detection integrator corresponds to the number of 
pulses returned from a target as the radar antenna scans 
through its beamwidth [Ref. 3]. For most practical radars, 
this number is sufficiently large so that the probability 
density function of the voltage at the output of the linear 
post-detection integrator shown in Figure 1 tends towards a 
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normal distribution which can be asymptotically expanded using 
Edgeworth's groupings of the Gram-Charlier series expansion 
[Ref. 9, Ref. 10]: 



f z {z) “ §{z') 



(29) 



<J>' 3) (z') 

6 



+ ^ 2 . *'4) u / } + 4>f6j (z /) 

24 ^ 72 Y 



_£L 4>(s) (z /) _ ik£k 4>n) (z /) 

120 144 



1296 



i — 4> (9) (z ; ) 



1 






where the terms on each line within the brackets are of the 
same order of magnitude [Ref. 11]; and 



z = integrator output voltage normalized by the rms 
noise voltage o at the input to the detector. 



az 



NK v.i 



V / ^ 



(30) 



$ (z) 



i r z : 
exp r — 



VTtT 



(31) 
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4> (m) (z) = (-l) m 4>(z) H m (z) 



( 32 ) 



H m ( z ) = mth order Hermite polynomial in z 




(32) 



The symbol K yr denotes the rth cumulant of the 
probability density function of the video filter output y 
which takes into account the additional noise-only samples 
averaged along with each signal-plus-noise sample. 

Studies by Cramer (Ref. 11] have validated the asymptotic 
behaviour of Edgeworth's series and determined that in 
general, the remainder term of the expansion is of the same 
order as the first term neglected. Hence, the output 
probability density function f z (z) can be expanded to the 
desired order of accuracy once the cumulants K y<r are known. 

The cumulants K y>r can be determined from Eqn (10) with X' 

y 

equal to — , i.e., 



K 



(2a 2 ) 1 (r-1) ! 



X 



[\fl + 2y 2 + 1] r - (Vl + 2 y 2 - 1] r 



(33) 




[ (1 +2y 2 ) 2 + l] r - [ (1 +2y 2 ) 



1 



_1 

2 




1 



1 



\ [ (1 + 2y 2 ) 2 + l] r + [ ( 1 + 2y 2 ) 2 - 1] r 
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Once the probability density function f 2 (z) is determined, 
the cumulative distribution function F z (z T ) can be found by 
integrating f z (z T ) in Eqn (30) from to z T . Hence, 



F z ( z T ) - [ P{z') 



(34) 



- 4> (2) (z') 

6 



+ — 1 (z ; ) + < ~ 1 ' <J> (5) (z f ) 



24 



72 



— <j> (4) ( z') - 



120 



£lSl 4>f6) (z /) _ _^_4>ca) (z /) 
144 1296 



where 



z / . qz r - ^y,l 






(35) 



•y.2 



P(^) “ /*' -p= exp [ 

J - c ° v/2tT l 2 



Z_ 2 1 
2 



j dz 



0.5 erf l 



\ 1 



z_ 

s[2 J 



(36) 



The probability of false alarm P £a and the probability of 
detection P d for a normalized threshold voltage of z T can then 
be determined by taking the complement of Eqn (34), i.e.. 



^Z, 0 ( Z T* 


(37) 


1 ^ 


(38) 
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where F ZiQ (z t ) denotes the cumulative distribution for the 
hypothesis that the target is absent while F ZiZ (z t ) denotes the 
cumulative distribution for the hypothesis that the target is 
present. 

The receiver operating characteristics is then obtained by 
plotting the probability of detection versus the probability 
of false alarm with X, y , and N as parameters. 

It is obvious from the above expressions that the process 
involved in the computation of the output probability density 
function f(z) and the receiver operating characteristics is 
numerically tedious. However, the steps involved are straight 
forward and can be easily programmed for execution in a 
digital computer. As part of this thesis, routines in PC- 
MATLAB have been written to compute the cumulants, the 
probability density functions, and the associated receiver 
operating characteristics. These routines are listed in 
Appendix B. 

Using these routines, the effect of insufficient video 
bandwdith on the the detection statistics of the Gaussian 
square-law system can be investigated for various values of y , 
input signal-to-noise ratio, and N. This is illustrated in 
Figure 3 to Figure 5. 

Figure 3 compares the output probability density functions 
obtained for a finite video bandwidth equal to a quarter of 
the IF bandwidth (Y“2) against the corresponding density 
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0.3 




Figure 3. Gaussian square law system output probability 
density function. 

functions when the video bandwidth is infinite ( y - 0 ) . The 
output probability density functions corresponding to the 
absence and presence of a target are respectively denoted by 
f z ,o(z) and f Zfl (z). For this illustration, the number of video 
samples integrated, N , is chosen to be nine and the signal-to- 
noise ratio in the IF range bin is assumed to be 3 dB when 
there is a target return. The choice of these parameters is 
in part to facilitate comparison with Marcum's plots [Fig. 4, 
Ref. 2] which assumed infinite video bandwidth. Accordingly, 
the plots in Figure 3 for y-0 are found to be in agreement 
with Marcum's results which serves as a validation of the 
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method and routines developed to compute the probability 
density functions. Graphically, the effect of the finite 
video bandwidth is observed to be two-fold. Firstly, the 
variance of f Zf0 (z) is reduced vis-a-vis the density function 
for infinite video bandwidth. Secondly, the variance of 
f Zrl (z) is also reduced but this is accompanied by a leftward 
shift of the mode. The net result is that for the same 
probability of false alarm, the detection performance has 
degraded considerably. This is the consequence of the 
integration of additional noise-only samples induced by the 
finite bandwidth of the video amplifier. 

Figure 4 shows the receiver operating characteristics for 
Y - 0 and y - 1 . The number of video samples integrated N is 
equal to 64 . For a given signal-to-noise ratio and 
probability of false alarm, it can be seen that the 
probability of detection achievable for the system with the 
finite video bandwidth (y=l) is poorer than the system with 
the theoretical reference of infinite video bandwidth. 
Conversely, the system with the finite video bandwidth will 
require a higher signal-to-noise ratio to have the same 
performance as the system with infinite video bandwidth. This 
is illustrated in Figure 5. For a probability of false alarm 
and probability of detection equal to 10' 7 and 0.95 
respectively, a system with infinite video bandwidth will 
require an input signal-to-noise ratio of 0 dB. To achieve 
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Probability of Detection Probability of Detection 




Figure 4. Gaussian square -law system operating 
characteristics 




Probability of False- Alarm 

Figure 5. Gaussian square-law system operating 
characteristics. 
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the same performance, the system with a video bandwidth that 
is half the IF bandwidth will require an input signal-to-noise 
ratio of 0.61 dB. From Eqn (1), this implies a collapsing 
loss of 0.61 dB that is attributable to the 'insufficient' 
bandwidth of the video amplifier. 
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III. 



COLLAPSING RATIO 



The statistical analysis in the preceding section provides 
the theoretical background and insight concerning bandwidth 
effects of a Gaussian-shaped video filter on the performance 
of the square-law detector. The analysis has also shown how 
sets of receiver operating characteristics such as those shown 
in Figure 4 can be generated using a digital computer to 
determine the input signal-to-noise required to attain a 
desired probability of detection for a specified probability 
of false alarm and a given IF bandwidth to video bandwidth 
ratio. While such an approach is feasible, it is more 
convenient in practice to follow Marcum's approach, using the 
collapsing ratio p defined in Eqn (2) to determine the 
required input signal-to-noise ratio. This is because the use 
of the collapsing ratio will obviate the need to generate sets 
of receiver operating characteristics for each value of y . 
As explained in Section I , only the receiver operating 
characteristics for an infinite video bandwidth (y-O) are 
necessary once the collapsing ratio is known. Furthermore, 
there are usually other sources of collapsing loss apart from 
insufficient video bandwidth in practical systems. In such 
cases, the formulation of the collapsing ratio induced by each 
source will facilitate computation of the overall collapsing 
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loss. A general treatment of the problem of multiple sources 
of collapsing loss can be found in Meyer [Ref. 4] and Barton 
[Ref 5]. It is however beyond the scope of this report. In 
this section, the aim is to formulate the collapsing ratio 
associated with insufficient video bandwidth for the Gaussian 
square -law system and to compare the result with Barton's 
formula for the collapsing ratio in Eqn (5). 

From Marcum's definition of the collapsing ratio, 

„ M + N 

P N (39) 

where M is the number of additional noise-only samples 
integrated along with the N signal-plus-noise samples. 

From the analysis of Section II, the number of additional 
noise-only samples integrated for every signal-plus-noise 
sample integrated in the Gaussian square law system is equal 
to N v -1 where N v is given by Eqn (26). It follows that if the 
number of signal-plus-noise samples integrated is equal to N, 
then 



M = N x (N v ~ 1) 



(40) 



and 



N(N-l) + N 



N„ 



P G - 



N 



(41) 



Hence, the collapsing ratio caused by insufficient video 
bandwidth in the Gaussian square-law system is equal to the 
equivalent number of independent IF samples N v averaged by the 
video filter. 

The subscript 'G' has been used to denote that the 
collapsing ratio p G given in Eqn (41) is caused by insuffi- 
cient video bandwidth in a Gaussian square-law system. It is 
also used to distinguish from Barton's formula for the 
collapsing ratio: 



Pa 



Bif + 2B V 

2B„ 



1 + y 



(42) 



For convenience, p a and p B shall respectively be referred 

as the Gaussian system formula and Barton's formula in the 
context of collapsing ratios caused by video bandwidth. As 
pointed out in Section I, Barton's formula models the actual 
filters with ideal rectangular shaped filters having band- 
widths equal to the noise-equivalent bandwidths of the actual 
filters. A comparison of the two formulae will provide a 
relative measure of the accuracy of Barton's formula in 
modelling systems that approximate the Gaussian square-law 
system. 

Figure 6 plots the collapsing ratio formulae p G and p B as 
a function of the ratio y of IF bandwidth to twice the video 

bandwidth. The plot obtained by applying Barton's formula is 
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Ratio of IF bandwidth to twice video bandwidth 

Figure 6. Collapsing ratio due to insufficient video 
bandwidth. 



indicated in the figure. As expected from Eqn (42), it is a 
linear plot and is independent of input signal-to-noise ratio. 
All the other curves are obtained by applying the Gaussian 
system formula for input signal-to-noise ratios of 0, 0.1, 1, 
and 2. It is is seen that p G is a function of the input 

signal-to-noise ratio X . For the same y / the collapsing 
ratio given by p G decreases as the signal-to-noise ratio 

increases. This dependency complicates the method of using 
the collapsing ratio to determine the required signal-to-noise 
ratio to attain a specified probability of detection and 
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probability of false alarm since the collapsing ratio is 
itself a function of the signal-to-noise ratio. For very 
small signal-to-noise ratios (X<1), the collapsing ratio p G 

can be approximated as 

Pc = Jl + 2y 2 (43) 

which is independent of the input signal-to-noise ratio. 

The collapsing ratio given by Eqn (43) can be considered 
as the worst case situation for any given y in the Gaussian 
square-law system. Taking p G -\/l + 2Y 2 as the reference, it is 
seen that Barton's formula overestimates the collapsing ratio 
for y less than two. In this region, the largest relative 
deviation is about 22.5% in the vicinity of y = 0-5. For y>2, 
Barton's formula underestimates the collapsing ratio. For 
large y, p B / p G approaches 1/^2” so that Barton's formula 
underestimates the collapsing ratio by a factor of about \[2 
for very large y . 

The collapsing loss associated with each collapsing ratio 
can be determined by finding the incremental increase in 
signal-to-noise ratio required with respect to the signal-to- 
noise ratio required in the ideal case of infinite video 
bandwidth for the same detection performance. Marcum has 
found that the collapsing loss is a function of the number of 
signal-plus-noise samples integrated but varies only slightly 
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with the probability of detection P a and probability of false 
alarm P ta [Ref. 2]. His results for the collapsing loss as 
a function of the collapsing ratio y and with N as a parameter 
is shown in Figure 7 for P d -0.5 and P fa -10" 10 . 




Coilopsing rotio p 



Figure 7. Collapsing loss vs. collapsing ratio for 
P d -0.5 and P fa - 10' 10 . (From Ref. 4) 



Using Figure 7 and the respective equations for p B and p G , 
the collapsing losses predicted by the application of Barton's 

formula are compared against those deduced using the Gaussian 
system formula for selected values of y in Table I and Table 
II. The collapsing loss associated with the use of Barton's 
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formula is denoted in the table by L C/b while that associated 
with the Gaussian system formula is denoted by L C/g . Using L Cig 
as the reference, the difference between the two collapsing 
loss for a given y is denoted by aL c . In Table I, the input 

signal-to-noise ratio X is assumed to be much less than one so 
that the collapsing ratio p B can be taken to be equal to 
v^l + 2y 2 . A signal-to-noise ratio equal to one is assumed in 

Table II to provide a relative bound on the collapsing loss 
for signal-to-noise ratios between zero and one. 

The comparison shows that for small signal-to-noise 
ratios, the collapsing loss predicted by Barton's collapsing 
ratio is within 0.3 dB of that deduced by the Gaussian system 
formula. If the signal-to-noise ratio is near or equal to 
one. Barton's formula overestimates the collapsing loss by 
about 0.4 dB for video bandwidths around half the IF 
bandwidth. 

Typically, the number of pulses integrated in a search 
radar may range from thirty to forty. In such cases, the 
required single-pulse input signal-to-noise ratio is close to 
unity for reliable detection performance [Ref. 5]. From 
Table II, it follows that the collapsing loss will be 
overestimated by about 0.4 dB if Barton's formula is used in 
a Gaussian square-law system with a video bandwidth about half 
the IF bandwidth. 
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TABLE I. COMPARISON OF COLLAPSING LOSS. 
BARTON'S FORMULA VS GAUSSIAN SYSTEM FORMULA. 



Input SNR 


Collapsing loss, dB 

C 




X « 1 


N • 10 


N - 1000 


Y 


Pm 


Po 


L c ,b 


L c.g 


AL c 


L c,b 


L 

c,g 


AL c 


1.0 


2.0 


1 .7 


0.7 


0.5 


♦0.2 


i . 1 


0.8 


♦0.3 


3.0 


4.0 


I 4.4 

| 


1.6 


1.8 


-0.2 


2.5 


2.8 


-0.3 


5.0 


6.0 


7.1 


2.2 


2.4 


-0.2 


3.3 


3.6 


-0.3 



TABLE II. COMPARISON OF COLLAPSING LOSS. 
BARTON'S FORMULA VS GAUSSIAN SYSTEM FORMULA. 



Input SNR 
X - 1 


Collapsing loss, L c dB 


N - 10 


N • 1000 


Y 


1 

Pm ! Po 

\ 

1 


L c,b 


L 

c,g 


AL c 


L c,b 


L 

c,g 


AL c 


1.0 


2.0 


1.5 


0.7 


0.3 


♦0.4 


1.1 


0.6 


0.5 


2.0 


3.0 


2.6 


1.2 


1.0 


♦0.2 


1.8 

i 


1.6 


i 

♦0.2 


5.0 


1 

6.0 


6.5 


2.2 

i 


2.3 


-0.1 


3.3 


3.4 


-0.1 



32 




IV. 



CONCLUSIONS 



The effects of Gaussian-shaped IF and video filters on the 
performance of a square-law detector in pulsed radar 
applications employing post-detection integration have been 
analyzed in this paper. 

Specifically, a formula for determining the number of 
additional noise-only samples integrated due to the finite 
bandwidth of the video filter has been derived. Emerson's 
formulation for the cumulants of the voltage probability 
density function at the output of the video filter has been 
used in this derivation. 

The resultant cumulants of the probability density 
function at the output of the post-detection integrator are 
obtained. Based on these cumulants and Edgeworth's asymptotic 
expansion of the Gram-Charlier series to represent the density 
functions, MATLAB routines have been written to plot the 
receiver operating characteristics for various ratios of IF 
bandwidth to video bandwidth. The degradation in performance 
or collapsing loss can be observed and quantified from these 
operating characteristics as the video bandwidth is reduced 
relative to the IF video bandwidth. 

Taking into account the Gaussian shape of the filters, 
the collapsing ratio is found to be equal to the equivalent 
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number N v of independent IF samples averaged by the video 
filter. This is given by Eqn (25) which is a quadratic 
equation in N v with signal-to-noise ratio and IF bandwidth to 
video bandwidth ratio as parameters. The collapsing ratios 
given by Eqn (25) and the associated collapsing loss are 
compared against those predicted by Barton's formula in Table 
I and II. The comparison shows that for typical signal-to- 
noise ratios close to 0 dB and video bandwidths about half the 
IF bandwidth. Barton's formula for the collapsing ratio 
results in an overestimation of the collapsing loss by about 
0.4 dB. 
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APPENDIX A 



MATHEMATICAL DERIVATION OF THE CUMULANTS 
This appendix follows the method of Emerson [Ref. 7] in 
deriving a closed form expression for the cumulants of the 
probability density function f Y (y) at the output of the 
Gaussian video filter shown in Figure A-l. 




Figure A-l. Block diagram of the receiver. 



The input to the system is a Gaussian random process given 



by : 



r 1 (t) = s ( t) + n ( t) 



(A-l) 



where 



sit) = s/2 A cos a > 0 t 



(A-2) 



and n(t) is additive white Gaussian noise ( AWGN ) with a power 
spectral density equal to tj /2 . 
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The impulse response functions of the Gaussian-shaped IF 



filter and video filters are respectively defined by: 




and 




( A- 4 ) 



where p and v are respectively the root-mean-square (rms) 
bandwidths of the IF and video filters. 

The starting point in this derivation is the following 
general expression of the cumulants K r obtained by Emerson for 
an arbitrary input signal contaminated by AWGN into a square- 
law detector with IF and video filters of arbitrary frequency 
response characteristics: 



“ (Tl)r (r ' 2 1) 1 E + (n> r- 1 r! Y, ^iSi(t ) 2 




( A-5 ) 



where 




(A-6) 
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X and ^(x) are respectively the eigenvalues and eigen- 



functions defined by the integral equation: 



X i $ i (x) = J g(x,y) ^(y) dy 



( A— 7 ) 



The function g(u,v) is called the system kernel and it is 
completely defined by the impulse response functions of the 
filters. In particular, it is defined by the following 
equation : 



which can be expanded into the uniformly convergent bilinear 
Mercer series, i.e., 



To determine the cumulants for the Gaussian system shown 
in Figure A-l, the method of solution is to evaluate the 
system kernel using Eqn (A-8) and expressing it in a form that 
matches the Mercer series. The eigenvalues and eigenfunctions 
that satisfy the integral equation of Eqn (A-7) can then be 
identified. The cumulants can then be determined by using 
these eigenvalues and eigenfunctions to evaluate Eqn (A-5). 




(A-8) 



g(u,v) = ^ X i 4 > i (u)<\> i (v) 



( A-9 ) 



i-o 
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Thus using Eqn (A-8), the system kernel for the Gaussian 



system shown in Figure A-l is given by 



g{u,v) - f 2y/2it P exp [ - (2?ip) 2 — — ] cos ( 0 o (u-z) 

- co 

y 2 

x y2jT v exp [ - (27tv) 2 — ] 

x 2^/2n P exp [ - (27cP) 2 ]coso> 0 (v- 

x dz 



Now, 

coso > 0 (u-z) coso > 0 (v-z) 
= -i-COS 0) o ( u- v) + 




Substituting Eqn (A-ll) into Eqn (A-10), 
=» g{u,v) = 47tP 2 /2lf [g x + g 2 ] 



z) 

( A-10 ) 
(A-ll) 



(A- 12 ) 
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where 



r \ 

g 1 - cost* > c (u-v) | exp ^ - 



(2m P ) 



[ ( u-z ) 2 + ( v-z ) 2 ] 



S^dl Z 2 ) ] dz 
2 J 



(A-13) 



^2 



f cos w 0 ( u+ v- 2 z) exp I - ( [ (u-z) 2 + (v-z) 2 ] 

J ^ 



+ (2nv)z z 2 ) ! dz 
7 J 



(A-14) 



The integral in Eqn (A-13) can be evaluated by completing 
the square in the exponent and casting the integrand into the 
form of the standard integral identity: [Ref. 12] 



f exp [-a 2 (x-m ) 2 ] dx - -^L 

J 3. 



(A-15) 



Thus completing the square, 



9 1 



COS <J> 0 ( u-v) 



exp |_ 



(2mp) 2 

2 



1+ V l (u 2 + v*)- 2£uv } 
l+2y 2 l+2y 2 



X 



J exp { 



^^( 1 + 2 y 2 ) [z l!_(u+v)] 2 ) dz 

Y 2 l+2y 2 



(A-16) 
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where 



Y 



- 1 - 



v 



B 



if 



2 B„ 



(A-17) 



The symbols JB 1£ and B v in the above equation are 
respectively the noise equivalent bandwidths of the Gaussian 
IF and video filters. 

04*2 ft 2 

Makinq use of Eqn (A-15) with a * — (l+2y 2 ) , 

Y 2 



9 1 



COS 0) o ( u-v) 



exp l 



(2nP) 2 

4 



x 



1 

_i 

i/2 t t v (l + 2y 2 ) 2 



{ — ( U 2 + V 2 ) 

1 + 2y 2 



2yhf_v j 
l+2y 2 ' 

( A- 18 ) 



Now, consider the second term g 2 of the system kernel 
given in Eqn (A-12). By completing the square in the exponent 
as in the evaluation of g x , it can be similarly shown that 



* - exp[-J^( 



l+Y 

1 + 2y 2 



( u 2 + v 2 ) - 






l + 2y'< 



( A— 19 ) 



oo 

x J cos co 0 ( u+ v-2z) exp [ - a 2 (z-m ) 2 ] dz 

— OO 



where 



a = (i+2y 2 ) 2 

Y 



7n - — ^ ( u+v) 

l + 2y 2 



( A- 20 ) 



( A- 21 ) 
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Now, 



j cosco 0 (u+v-2z) exp [ -a 2 (z-m ) 2 ] dz 



(A 



J cos co 0 ( u+v) cos 2co 0 z exp [ -a 2 (z-m ) 2 ] dz 



+ J sinco D ( u+ v) sin2io 0 z exp [ -a 2 (z-m ) 2 ] dz 



Let x = (z-m) . Then, 



J cos<o 0 (u+v) cos2to 0 z exp [ -a 2 (z-m ) 2 ] dz 



(A 



cos co 0 ( u+v) 003 2 ( 0 ^ J cos 2 ( 0 o x exp [ -a 2 x 2 ] dx 

co 

cos co 0 ( u+v) 3102(0^ j sin2oo 0 x exp [ -a 2 x 2 ] dx 



From Eqn (679) of Ref. 12, 



J cos2o> 0 x exp [ -a 2 x 2 ] dx - exp [ 






a 2 J 



(A 



J sin2o> 0 x exp [ -a 2 x 2 ] dx - 0 



(A 



- 22 ) 



-23) 



-24) 

-25) 
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Hence , 



J cos(i> 0 (u+v) cos2(i) a z exp[-a 2 (z-iii) 2 ] dz 



(A-26) 



fr r o>i 

cos 0 > o (u+v) cos2 wjn -^-expj ^- — - 



The second integral in Eqn (A-22) can be similarly 
evaluated to give the following result: 



J si no) a (u+v) sin2 <o c z exp [ -a 2 (z-m) 2 ] dz 



(A-27) 



[r r wi 

sinto (u+v) sin2tojn -*■ — exp i 

a 1 a 2 



Using the results of Eqn (A-22), Eqn (A-26), and 
Eqn (A-27) to evaluate Eqn (A-19), it can be shown that 



g 2 . exp! ( JLllL. (u2+v 2).^1!£Z )j 

L 4 l+2y 2 l+2y 2 1 



( A-28 ) 



exp l 






27t 2 p 2 (l+2y 2 ) 



cos { o>.(u+v) — } 

l + 2y 2 

_i 

y/2n v (l+2y 2 ) 2 



Comparing Eqn (A-28) with Eqn (A-18), it is obvious that 
if to 0 > 2 tiP, then g 2 will be small and negligible compared to 
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Hence, neglecting the contribution from g 2/ the system 



9 1 - 

kernel given in Eqn (A-12) becomes 



g(u, v) 



47rP 2 COS0) o (u-v) 
y/l + 2y z 



( A-29 ) 



x 



exp y 



(2np) 2 

2 



{ — ( u 2 + v 2 ) 
l + 2y 2 



2 y 2 u v j 
l + 2y 2 J 



Using the 


following 


substitutions , 




u = 


i 

( l + 2y 2 ) 4 


v 


(A— 30 ) 


2 7T P 






v = 


_1 

(l+2y 2 ) 4 


y 


( A- 31 ) 


2 7t P 




t = 


v / l+2y 2 -1 

VT+2y 2 +1 




(A-32) 



Eqn (29) can be transformed into the following form suitable 
for the application of Mehler's formula: 

g(u,v) - [ cosu 0 u coso) 0 v + sino) c u sino> 0 v ] (A-33) 

4tiP 2 [ _ x 2 +y 2 1 f _ t 2 (x 2 +y 2 ) - 2 txy ' 

1 J 1 J 

Using Mehler's formula, 



exp y 



\._ tHx^y2) -2txy H,U)H,(y) 






i-o 



i! 2 J 



( A-34 ) 
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where H L (x) is the Hermite polynomial defined by 



H i (x) = ( -1 ) i exp [x 2 ] — , exp[-x 2 ] 

dx J 

From Eqn (A- 30) and Eqn (A-31), 
x = 2tiPu 

J. 

(l+2y 2 ) 4 

y = 

_1 

( l+2y 2 ) 4 

Substituting the above equations into Eqn (A-33), 

g(u,v) = g c (u,v) + g s (u,v) (A-38) 



( A-36 ) 
( A- 3 7 ) 



where 



g c (u, v ) 



2Py^T j- y/i+2y 2 - 1 
1 + v / l+2y 2 v/l+2y 2 +1 



x 



x 



X 



\/4tTP exp [ - ( — j cos w 0 u 
v / l+2y 2 2 

_i _i 

[ 2 1 i ! v/tT (l+2y 2 ) 4 ] 2 

v/4t:P exp [ - J cosw 0 v 

L y/l^F 2 

_1 _1 

[ 2 i i! v/tT (l + 2y 2 ) 4 ] 2 

E rr \ 2tc 6 u 1 „ f 2 tc 6 v 

[ li — T J ^ l E T 

( l + 2y 2 ) 4 ( l + 2y 2 ) 4 



( A-39 ) 



44 



( A-40 ) 



g s (u, v) 



2py/rc~ f y/l+2y 2 -1 ] i 



1 +v/l + 2y 2 v ^1 + 2 y T + 1 



/t — it- f (2pe6) 2 u 2 1 

v/ 4W1T exp I - v & t — i — — j coso> 0 ii 

y/l+2y 2 2 



1 1 



[ 2 i i ! v/i - (l+2y 2 ) 4 ] 2 



7i |i exp [ - ( 27t P )2 - J sinw 0 v 

yl + 2y 2 z 

_1 _1 

[ 2 i i ! v/tT ( l+2y 2 ) 4 ] 2 



E M 



2 uP u 1 „ r 2 uP v 



i-0 



r J^L 



(l + 2y 2 ) 



2 \ 4 



(l + 2y 2 ) 



2\ 4 



By comparing Eqn (A-38) against the Mercer series 
expansion of the system kernel 



g c {u,v) - £ 4> C/i (v) 



( A- 41 ) 



i-0 



the eigenvalues and eigenfunctions can be deduced to be 



C , 1 



2py/F r y/l4-2y 2 -1 ^ 

i+ Vi+2y 2 Vi+2y 2 +i 



( A-42 ) 



v/iirpr exp [ - 






(2llP ) 2 U 2 1 rr I” 27tP U 

2 J 1 L a J 

v 1+2Y (l+2y 2 ) 4 



cosw 0 u 



_i _i 

[ 2 2 i! </tz (l+2y 2 ) 4 ] 2 



( A- 4 3 ) 
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Similarly, it can be shown that 



K. 



4> s 



Hence 

oo 

E 

i-0 



Now, 

Hence 



k 



c , i 



( A- 4 4 ) 



tTP" exp y 



f ( 2 TC P ) 2 u 2 
2 



ri (u) 



v / T+2y 1 



1 w ^ 



2nftu 



( 1 + 2 Y 2 ) 4 



sina) 0 u 



[ 2 i i ! ,/F (l + 2y 2 ) 4 ] 2 



, from Eqn (A-42) and Eqn (A-44), 



( A-45) 



(^c,i) r - E (X s,i) r (A~46) 

i-0 

(2Py /F) r |- \/l + 2Y 2 ~ 1 ] ri 

(i +yi+2 Y 2 ) r i_ ° yi+2 Y 2 + 1 

CO 

Y' x J = — - — for I x I < l . 

4 -^ 1 - V 



oo 



E 

i-0 




( 2 Py/F) r 

(V1 + 2Y 2 +D r - (V1 + 2Y 2 -1) r 



GO 



E < l i )r 



E * E <*..*>' 



i-0 



i-0 



( A-47 ) 



2(2p/F) r 

(v/l + 2Y 2 +l) r - (yi + 2 Y 2 -l) r 



(A-48) 
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From Eqn (A-12) and Eqn (A-45), 



s ci ( t) = f2A J cos u> 0 ( t-u) 4> Cji (u) 



du 



(4rcP) 2 



_i _i 

[ 2 i i ! s/W (l + 2y 2 ) 4 ] 2 



(A-49) 



/ LZ ^ 

exp(-a 2 — ) tf^au) cos co 0 u cos u) 0 ( t-u) 
2 



du 



where 



2kP 

_i 

(l + 2y 2 ) 4 



( A-50 ) 



1 1 

COS 0) o u COS 0) o ( t-u) = — COS( 0 o t+ — COS (20) o u+ t) 

2 2 

( A- 51 ) 

For co 0 > 2 tcP , the residue from the high frequency 

component in Eqn (A-51) can be neglected in the integration of 
Eqn (A-49). Hence, 



y/2A(47Cp) 2 CQSC0 o t (A- 5 2 ) 

11 2 

[ 2 i i\ fit (l + 2y 2 ) 4 ] 2 

oo 

x J exp ( -a 2 —- ) H i (au) du 

- oo 
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The following integral identities derived in Ref. 13 are 



useful in evaluating (A-52): 



J exp(-x 2 ) H 2n+1 {j2x) dx - 0 



( A- 53) 



/ exp ( -x 2 ) H z Jy/ 2x) dx = (2nUJW 



( A-54 ) 



Let a 2 — = x 2 

2 



(A-55) 



Then au = JTx and du - dx . 

a 

Using the above substitution, Eqn (A-53) and Eqn (A-54) 



J exp(-a 2 -^-) H i (au) du 
= ^ - ^ 2 p Y ^ — j exp ( -x 2 ) H i ( s/2x) dx 



r v/2tT ( l+2y 2 ) 4 i ! 

2n P , for i even 

i 2 

( A-56 ) 

l 0 for i odd 
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Hence, using the result of Eqn (A-56) to evaluate Eqn (A-52) 



( 

2 P/rT 



(l + 2y 2 ) 



(i!) 2 



COS (i) 0 t 



s 



,i<t> 




0 



for i odd. 



for i even. 



( A-57 ) 



A similar result is obtained for s Sti (t) except that the 
coso) 0 t term in Eqn (A-57) is replaced by a sinto 0 t term. 

Using Eqn (A-40) and Eqn (A-55), it can be shown that 



E u c ,i) r 4i(t) 



2A 2 ( 2P_y/i V ^ ( i+2y 2 ) 4 cos 2 0 ) o t 
( l+/l+2v r ) r 



j- y /, l + 2y 2 ~ 1 ] 2 ri (2 i) ! 

— \Zl + 2y 2 +1 2 2i (i ! ) 2 



Let [ -^ 1+2y2 -1 ] 2r - 
V r l+2v 2 +1 






( A-58 ) 
( A-59 ) 



From Ref . 12 , 



r i (2 i) ! 



1-0 



2 2i ( i ! ) 2 



1 + — x + 
2 



4 (4 + Dx 2 
2 2 

2 ! 



(1-X) 2 



( A- 60 ) 
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Hence, using this result to evaluate Eqn (A-58) gives 



E <*c.i> r 4i<t> 



_i 

2A 2 (2 Py/rT) r ~ 1 (l + 2y 2 ) 4 

_1 

[ (V1+2Y 2 +1) 2r - (Vl+2y 2 -1) 2r ] 2 



COS 2 0) o t 



( A- 61 ) 



Similarly, it can be shown that 



E 4i(t> 



j. 

2A 2 (2p/F ) r-1 (l + 2y 2 ) 4 

_i 

[ (Vl + 2y 2 +1) 2r - (V1 + 2Y 2 -1) 2r ] 2 



sin 2 0 ) o t 



Hence , 



( A-62 ) 



E ^i) r sUt) 

i-0 



E (X c,i) r + E ^s.i) r sli(t) 



_1 

2 A 2 (2Py/7T) r ~ x (1 + 2Y 2 ) 4 

_1 

[ (v/1 + 2y 2 +l) 2r - (V 1+2Y 2 -l) 2r ] 2 



( A-63 ) 
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From Eqn (A-5) and the results of Eqn (A-48) and Eqn (A-63), 



K z = (r\)* . < r f±l £ (X ,) 1 + ( T| ) r_1 •t ! £ n z ) r s?(t) 

Z 2-0 2-0 



( 2r| P s/tT ) r (r-1) ! 
(Vl+2y 2 +D r - (\/l + 2Y 2 -D r 



r 

L 



l+r- 






(y/l + 2Y 2 ->-l) r - (\/1 + 2 y 2 -l) r 



N 



[ (v^l + 2 Y 2 +1 ) r - (Vl + 2Y 2 -l) r ] 



JL 

2 



X 



(1+2y 2 ) 




( A— 64 ) 



Now, it can be shown that the IF noise bandwidth B if is 
equal to py/iT. It follows that the total average power into 
the detector is given by 

o 2 - T|P/n~ (A— 65 ) 



Hence, 



K r 



(2o 2 ) f (r-1) ! 

(Vl + 2Y 2 +l) r - (v/l + 2Y 2 -l) r 



( A- 6 6 ) 



r 

i 



i 



r jii (y^2^ + D r - (y^2^-D f (1+2Y 2)i] 
° 2 \ ( v / l + 2Y 2 +l) r + (Vl+2Y 2 -D r J 
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This equation which expresses the cumulants K r in terms 



of the IF bandwidth to video bandwidth ratio y 
signal-to-noise ratio A 2 / a 2 concludes the 
derivation of this appendix. 



and the input 
mathematical 
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APPENDIX B 



MATLAB PROGRAMS FOR PLOTTING 
THE RECEIVER OPERATING CHARACTERISTICS 

This appendix is a listing of the Matlab programs written 
as part of this thesis to plot the probability density 
functions and operating characteristics of the Gaussian square 
law system described in Section II. The programs are written 
as MATLAB M-Files. The following are respectively the 
filenames of the programs for plotting the probability density 
functions and operating characteristics: 

1. SQ2PDF .M: This plots the probability density functions 

based on Eqn (29) which is an Edgeworth's asymptotic 
expansion of the Gram-Charlier series. 

2. SQ2ROC .M: This plots the receiver operating characteris- 

tics based on Eqn (37) which is derived by integrating 
Eqn (29) from the threshold voltage to infinity. 

These are script files which are run by invoking their 
filenames in the MATLAB operating environment. Both programs 
will prompt the user to enter the necessary input parameters 
prior to execution. In addition five function type M-Files 
are specifically created to facilitate programming and 
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execution of these script files. The filenames of these five 
functions and their purpose are: 

1. GCUMU.M : This function computes the cumulants of the 

video filter output probability density funtion based on 
Eqn (10). 

2. GDFZ.M : This function computes the probability density 

for given normalized values of the Gaussian square-law 
system output based on Eqn (29). 

3. UPFZ.M : This function computes the probability of false 

alarm or probability of detection of the Gaussian square- 
law system for a specified threshold voltaqe. The 
computation is based on Eqn (37) and Eqn (38). 

4. GCR.M : This function computes the collapsing ratio 

caused by the finite video bandwidth. It is based on Eqn 
(26). 

5. HMT.M : This function computes the mth order Hermite 

polynomial for values of m less than or equal to 12. 



The MATLAB source programs of the two script files and the 
five function files are provided in the remaining pages of 
this appendix. 
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1. SQ2PDF.M PROGRAM 



-6 

% This program computes and plots the pdf of a square law 
% detector with Gaussian IF and video filters. It will 
% prompt the user to enter the following parameters: 

% a. Range and resolution of the x-axis for the plot. 

% b. SNR = input signal-to-noise ratio. (Not in dB) 

% c. GAMMA = IF bandwidth to twice video bandwidth ratio. 

% d. N = Number of pulses to be integrated. 

% In addition to the pdf plot for the specified GAMMA, the 
% pdf plot for an infinite video bandwidth will automatically 
% be plotted for comparison purposes. 



disp('The output pdf will be plotted as a function of the' 
disp( 'normalized voltage, z.') 

disp( 'Enter range and resolution of z to be plotted') 
range = input ( ' [ zmin, zmax, zres] : ') 

N = input ( 'Number of pulses to be integrated = ') 

SIGMA = input (Input noise rms voltage = ') 

SNR =[00] 

GAMMA =[00] 

ORDER = 4 

SNR(2) = input ('Enter input SNR: ') 

GAMMA ( 2 ) = input ('Enter value of GAMMA: ') 

disp( 'Computing ...') 

KYO = [ ] 

KYI = [ ] 

Z = range ( 1 ) : range ( 3 ) : range ( 2 ) 

GAMMAO = GAMMA ( 1 ) 
for j = 1:2 
K = [ ] 

SNR j = SNR(j) 
for r = 1: ORDER 

K(r ) = gcumu (r, SIGMA, SNR j, GAMMAO) 
end 

if j == 1 

fZinfO = gdpz ( ORDER, N,Z, SIGMA, SNR j, GAMMAO) 
else 

fZinfl = gdpz ( ORDER, N,Z, SIGMA, SNR j, GAMMAO) 
end 
end 

GAMMA V = GAMMA ( 2 ) 

SNR j = 0.0 

for r = l:ORDER+2 

KYO ( r ) = gcumu (r, SIGMA, SNR j,GAMMAv) 
end 

f ZvO = gdpz ( ORDER, N,Z, SIGMA, SNR j,GAMMAv, KYO) 
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SNR j = SNR( 2 ) 

NPRIME = gcr ( SNR j , GAMMAv ) 

SNRjbar = SNR j/NPRIME 
for r = l:Order+2 

KYl(r) = gcumu( r , SIGMA , SNRjbar , GAMMAv ) 
end 

fZvl = gdpz ( ORDER, N,Z, SIGMA, SNR j, GAMMAv, KYI) 

plot ( Z , f Zinf 0 , '-' ,Z,fZvO, ' — ' , Z , f Zinf 1 , '-' ,Z,fZvl, ' — ' ) 

grid 

xlabel ( 'Normalized Output Voltage, z') 
ylabel ( ' Probability Density, f(z)') 

% ===================================================== 



2. SQ2ROC.M PROGRAM 



% ========================================================== 

% This program computes and plots the ROC data points for 
% the Gaussian square law system. It is programmed to provide 
% accuracy up to the 4th grouping of Edgeworth's series. To 
% obtain any other order of accuracy, the local variable 
% ORDER which is set to 4 and the function file gdfz.m must 
% be modified accordingly. 

% 

disp( 'Choices: ' ) 

disp('l. Computing Pfa and Pd given a normalized threshold.') 
disp('2. Pd vs Pfa curves.') 
disp( ' ') 

CHOICE = input ( 'Select one of the above: ') 

if CHOICE == 1 
ORDER = 4 
MORE = 1 
while MORE == 1 

N = input ( 'Number of samples integrated = ') 

SNR = input ('Input SNR = ') 

SIGMA = input ('Input rms noise voltage = ') 

GAMMA = input ( 'GAMMA = ') 

VT = input( 'Normalized threshold voltage, VT = ') 

disp( 'Computing ... ') 
for r = 1 : ORDER+2 

KYO(r) = gcumu(r , SIGMA, 0, GAMMA) 
end 

NPRIME = gcr (SNR, GAMMA) 

SNRjbar = SNR j/NPRIME 
for r = 1: ORDER+2 

KYl(r) = gcumu(r, SIGMA, SNRjbar , GAMMA) 
end 
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Pfa = upf z ( ORDER , N , VT , SIGMA , KYO , 1 ) 

Pd = upf z( ORDER, N,VT, SIGMA, KYI, 1) 

disp( 'Do you wish to compute Pd vs Pfa for another VT?') 
MORE = input ( 'Type 1 if YES, 0 if NO : ') 

end 
else 

N = input ( 'Number of samples integrated = ') 

GAMMA = input ('GAMMA = ') 

disp( 'Enter up to a maximum of 3 SNR values as elements of 

[]') 

SNR = input ( ' SNR [] = ') 

SIGMA = input( 'Input rms noise voltage = ') 

disp( 'Enter range & resolution of z to be used in 
computing' ) 

disp('the ROC.') 

range = input ( ' [zmin, zmax, zres] ' ) 

Z = range ( 1 ) : range ( 3 ) : range ( 2 ) 
for r = 1 : ORDER+2 

KYO(r) = gcumu (r, SIGMA, 0, GAMMA) 
end 

NPRIME = gcr( SNR, GAMMA) 

SNRjbar = SNR ./NPRIME 
for r = 1: ORDER+2 

KYl(r) = gcumu(r, SIGMA, SNRjbar , GAMMA) 
end 

Pf a ( j , : ) = upf z (ORDER, N,Z, SIGMA, KYO, 0) 

Pd ( j , : ) = upfz(ORDER,N,Z, SIGMA, KYI, 0) 
if length (SNR) == 1 
semi logx( Pfa ,Pd , '-' ) 
elseif length(SNR) == 2 

semilogx(Pf a( 1 , : ) , Pd ( 1 , : ) , Pf a ( 2 , : ) , Pd( 2 , : ) ) 
elseif length(SNR) == 3 

semi logx( Pfa ( 1 , : ) ,Pd( 1 , : ) , Pf a( 2 , : ) , Pd( 2 , : ) , Pf a( 3 , : ) , Pd( 3 , : ) ) 
end 
grid 

xlabel( 'Probability of False Alarm') 
ylabel( 'Probability of Detection') 
end 
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3. GCUMU.M FUNCTION 

function K = gcumu(r, SIGMA, SNR, GAMMA) 

% ========================================================== 

% This function computes the rth cumulant of the output pdf 
% of a square law detector with Gaussian IF and video filters. 
% GAMMA is the ratio of the IF bandwidth to twice the video 
% bandwidth. SIGMA is the input rms noise. 

% Reference: Thesis Eqn (10). 

% 

G12 = 1 + 2 *GAMMA A 2 
G12PR = (sqrt(G12) + l) A r 
G12MR = ( sqrt(G12 ) - l) A r 
BN1 = (2 * SIGMA A 2) A r 
RM1 = r - 1 
BN 2 = fac(RMl) 

BD1 = G12PR - G12MR 
B = ( BN1 * BN2 ) /BD1 
CN1 = BD1 

CD1 = G12PR + G12MR 
C = sqrt ( CN1/CD1 ) 

K = B * (1 + r*SNR*G12 A 0 .25) 

% ============== ========================================= 



4. GDFZ.M FUNCTION 

function fz = gdf z( ORDER, N,Z, SIGMA, SNR, GAMMA, K) 

% ========================================================== 

% This function computes the probability density for specified 
% values of the Gaussian square law system normalized output, 
Z. 

% K is an array comprising the cumulants of the pdf. 

% Reference: Thesis Eqn (29). 

% 

C = [ ] 

for j = 1: ORDER 
SUMn =0.0 
SUMd =0.0 
for m = 1:N 

SUMn = SUMn + K(m,j+2)/N 
SUMd = SUMd + K(m, j+2 )/N 
C(j) = SUMn/(N A ( j/2) * SUMd A ( ( j+2)/2) ) 
end 

fz = zeros ( 1 , length(Z) ) 

SUM = zeros (ORDER, length(Z) ) 

X = ( SIGMA*Z - N*K(1) )/sqrt(N*K(2) ) 

CONST = SIGMA/sqrt(N*K( 2 ) ) 
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PSIO = exp(-(Z . A 2 )/2 ) /sqrt ( 2*pi ) 
for j = 1: ORDER 

if j == 1 

SUM( 1 , : ) = CONST * PSIO .*( 1 + hmt(Z,3) * C(l)/6 ) 
elseif j == 2 

SUM( 2 , : ) = CONST * PSIO .*( hmt(Z,4) * C(2)/24 + hmt(Z,6) 

* (C( 1 ) ) A 2/72 ) 

elseif j == 3 

SUM( 3 , : ) = CONST * PSIO .*( hint ( Z , 5 ) *C ( 3 ) /12 + hmt(Z,7)* 
C(1)*C(2)/144 + hmt( Z,9)*(C(1) ) A 3/1296 ) 

elseif j == 4 

SUM( 4 , : ) = CONST * PSIO .*( hmt ( Z , 6 ) *C ( 4 ) /720 + hmt(Z,8)* 
( ( C ( 2 ) ) A 2/1152 + C( 1 ) *C( 3 ) /720 ) + hmt(Z,10)* 
( C ( 1 ) ) A 2 * C( 2) / 1 7 2 8 + 

hmt ( Z , 1 2 ) * ( (C(l) ) A 4/31104 ) 
end 

fz = fz + SUM( j , : ) 
end 



5. UPFZ.M FUNCTION 

function UPTE = upfz ( ORDER, N , Z , SIGMA, ERROR) 

% ========================================================== 

% This function computes the Pfa or Pd of a Gaussian square 
% law detector for a specified threshold normalized by the rms 

% noise voltage at the input to the detector. If the argument 
% ERROR is set to 1, the function will return a two element 
% row vector. The second element indicates the order of 
% accuracy of the Pfa or Pd value given by the first element. 

% If ERROR is set to 0, then the function will return the 

% the value of Pfa or Pd without any indication of the order 
% of accuracy. 

% Reference: Thesis Eqn (37) and Eqn (38). 

% 

C = [] 

for j = 1: ORDER 
SUMn = O.o 
SUMd =0.0 
for m = 1:N 

SUMn = SUMn + K(m,j+2)/N 

SUMd = SUMd + K(m,2)/N 

end 

C ( j ) = SUMn/(N A ( j/2 ) * SUMd A ( ( j+2)/2) ) 
end 

ZZ = (SIGMA*Z - N*K(1) )/sqrt(N*K(2) ) 

PSIO = exp(-ZZ . A 2 )/sqrt( 2*pi ) 

ZM = ZZ/sqrt ( 2 ) 

SUM1 = 0 . 5*erf ( -INF , ZM , "high 7 ) 
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TERM1 = hurt ( ZZ , 2 ) *C( 1 ) / 6 

TERM2 = hmt ( ZZ , 3 ) *C ( 2 ) /24 + hmt( ZZ , 5) *C( 1 ) A 2/72 

TERM 3 = hmt (ZZ,4)*C(3)/120 + hmt ( ZZ , 6 ) *C ( 1 ) *C ( 2 ) /144 + 

hmt(ZZ,8)*C(l) A 3/1296 

TERM4 = hmt( ZZ,5)*C(4)/120 + hmt ( ZZ , 7 ) *C ( 2 ) A 2/1152 + 

hmt (ZZ,7)*C(l)*C(3)/720 + hmt ( ZZ,9)*C(1) A 2*C(2)/1728 
hmt( ZZ,11)*C(1) A 4/31104 
SUM2 = PSIO . * ( TERM1 + TERM2 + TERM3 + TERM4 ) 
if ERROR == 1 

UPTE(l) = 1 - ( SUM1 - SUM2 ) 

UPTE ( 2 ) = PSIO . *TERM4 
else 

UPTE = 1 - ( SUM1 - SUM2 ) 
end 

% ========================================================== 



6. GCR.M FUNCTION 
function GRHO = gcr ( SNR , GAMMA ) 

% ========================================================== 

% This function computes the collapsing ratio caused by 
% insufficient video bandwidth in a Gaussian square law 
% system. GAMMA can be a vector but SNR must be single 
valued. 

% Note that GRHO is also the equivalent number of independent 
% samples averaged by the video filter. 

% Reference: Thesis Eqn (26). 

% 

B12 = 1 + 2*GAMMA . A 2 
Bll = 1 + GAMMA . A 2 

B = 2*SNR*sqrt ( B1 2 ./Bll) - sqrt(B12) 

C = -2*SNR*sqrt ( B12 ) 

GRHO = 0 . 5* ( -b + sqrt ( B . A 2 - 4*C) 

% ========================================================== 



7. HMT.M FUNCTION 
function H = hmt(Z,m) 

% ========================================================== 

% This function computes the mth order Hermite polynomial for 
% values of m less than or equal to 12. 

% 

if m == 0 
H = 1 

elseif m == 1 
H = Z 
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2 



elseif m == 

H = Z . A 2 - 1 
elseif m == 3 

H = Z . A 3 - 3*Z 
elseif m == 4 



H = Z . A 4 - 6*Z . A 2 + 3 
elseif m == 5 



H = Z . A 5 




10*Z 


. A 3 


+ 


15*Z 








elseif m == 


6 
















H = Z . A 6 


- 


15*Z 


. A 4 


+ 


45*Z . 


A 2 ■ 


- 15 




elseif m == 


7 
















H = Z . A 7 


- 


21*Z 


. A 5 


+ 


105*Z 


. A 3 


- 105*Z 




elseif m == 


8 
















nc 

ii 

N 

• 

> 

OD 


- 


28*Z 


. A 6 


+ 


210*Z 


. A 4 


- 420*Z . A 2 + 


105 


elseif m == 


9 
















H = Z . A 9 


- 


36*Z 


. A 7 


+ 


378*Z 


. A 5 


- 1260*Z . A 3 + 


945*Z 


elseif m == 


10 
















H = Z . A 10 


- 


45*Z 


. A 8 


+ 


630*Z 


. A 6 


- 3150*Z • A 4 + 


4725*Z 


- 945 


















elseif m == 


11 
















H = Z . A ll 


- 


55*Z 


. A 9 


+ 


990*Z 


. A 7 


- 6930*Z . A 5 + 


17325*Z 



- 10395*Z 
elseif m == 12 

H = Z . A 12 - 66*Z . A 10 + 1485*Z . A 8 - 13860*Z . A 6 
+ 51975*Z . A 4 - 62370*Z . A 2 + 10395 

end 

% ================================================= 
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